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Abstract: We investigate the cosmological aspects of Tree Level Gauge Mediation, a 
recently proposed mechanism in which the breaking of supersymmetry is communicated to 
the soft scalar masses by extra gauge interactions at the tree level. Embedding the mecha- 
nism in a Grand Unified Theory and requiring the observability of sfermion masses at the 
Large Hadron Collider, it follows that the Lightest Supersymmetric Particle is a gravitino 
with a mass of the order of 10 GeV. The analysis in the presence of R-parity shows that 
a typical Tree Level Gauge Mediation spectrum leads to an overabundance of the Dark 
Matter relic density and a tension with the constraints from Big Bang Nucleosynthesis. 
This suggests to relax the exact conservation of the R-parity. The underlying SO (10) 
Grand Unified Theory together with the bounds from proton decay provide a rationale for 
considering only bilinear R-parity violating operators. We finally analyze the cosmologi- 
cal implications of this setup by identifying the phenomenologically viable regions of the 
parameter space. 
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1 Introduction 

Low-energy supersymmetry (SUSY) provides an appealing framework for the origin of the 
Dark Matter (DM) component of the energy density of the Universe. In the Minimal 
Supersymmetric Standard Model (MSSM) only two particle states stand out for their 
phenomenological viability: the neutralino and the gravitino. 

In this work we investigate the case where the DM candidate is the gravitino with a 
mass in a specific range dictated by Tree Level Gauge Mediation (TGM) [1-3], a recently 
proposed mechanism in which the breaking of SUSY is communicated to the soft scalar 
masses by extra gauge interactions at the tree level. 

The mass of the gravitino is particularly sensitive to the mechanism of SUSY breaking 
and its mediation to the visible sector. In a large class of models SUSY is broken by the 
F-term of a chiral superfield {Z) = F9^ and the breaking is communicated to the Standard 
Model (SM) chiral superfields at the mediation scale M. The two most popular mechanisms 
of mediation of SUSY breaking are gravity and (loop) gauge mediation. 

In the former case the scale M coincides with the Planck scale, Mp. Then the soft 
terms and the gravitino mass are expected to be of the order of F/Mp, implying that the 
gravitino is not always the LSP and its mass is comparable with mgoft ~ 1 TeV. On the 
other hand, in theories based on loop gauge mediation the soft scale is 



a F 



(1.1) 
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depending on the two parameters F and M. In this case the gravitino is always the LSP 
with a mass ranging from 0{eV) to O(GeV). 

In this paper we investigate the cosmology of the alternative scenario given by TGM. 
In such a case the chiral superfield Z is a SM singlet with charge Xz under U{l)x, which 
is an extra abelian gauge symmetry broken at the scale Mx- The MSSM superfields Q 
carry a non trivial charge Xq under U{l)x- Integrating out the heavy vector superfield 
Vx related to the U{l)x factor, soft masses are induced through the F-term (Z) = FO"^ 




Figure 1: TGM super-graph generating the bilinear soft masses. 

and read 

(?^Q)trcc = g\XQXzj^ , (1.2) 

where gx is the gauge coupling relative to U{l)x- 

The gravitino mass is then related to mgoft by the simple relation 

Mx , . 

m3/2 ~ ^soitjj^ ■ (1-3) 

It is natural to embed the the enlarged Gsm <S5 U{1)x group into a rank-5 Grand Unified 
Theory (GUT) such as (5*0(10). Hence we assume that the SUSY breaking is communicated 
at the S'O(IO) — Gsm breaking scale Mq ~ 10^^ GeV. We note a few interesting features 
of this model: 

• the sfermions masses are flavour universal as in (loop) gauge mediation, thus solving 
the supersymmetric flavour problem. 

• from Eq. (1.2) we can see that the ratios among the soft masses {mQ)tvee depend just 
on the choice of the U{l)x charges, making the model predictive and testable at the 
LHC. 

• given Mx = Mq, Eq. (1.3) implies that the gravitino mass is of the order of 10 GeV^. 

The outline of the paper is the following: in the next section we present the most relevant 
features of TGM and set the framework for the cosmological analysis. In the presence of 
exact R-parity the outcome is that for a typical TGM spectrum the DM relic density is 
overabundant and the decay of he NLSP is in tension with the BBN constraints. On the 

^Such value for the gravitino mass is not typically accessible in standard (loop) gauge mediation. On 
the other hand LSP gravitinos of 0(10 GeV) are also achievable in scenarios like gaugino mediated SUSY 
breaking [4, 5]. 
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other hand a smah amount of R-Parity Violation (RPV) can easily restore the agreement 
with the BBN and, at the same time, reduces the amount of gravitinos produced by the 
NLSP decay. Hence we study how RPV can be obtained in the context of an 50(10) model 
of TGM and we show that the GUT structure of the theory motivates the restriction 
to bilinear RPV. We finally analyze the cosmology of such a scenario and identify the 
phenomenologically viable regions of the parameter space. 



2 5*0 (10) Tree Level Gauge Mediation 

As already mentioned in the introduction we are going to analyze the 50(10) TGM model 
presented in [1]. We provide here a brief overview of the main features of the low-energy 
TGM spectrum. 

The U{l)x gauge group responsible for the TGM mechanism is identified with the 
abelian factor external to SU{5) in the embedding SU{5) U{l)x ^ 50(10). After the 
one-step breaking 50(10) — > Gsai at the scale Mq, all the effects of TGM and the GUT 
physics are encoded in the MSSM boundary conditions at the GUT scale. 

The ratios among the tree level soft masses (?fig)tree depend only on the embedding of 
the MSSM chiral superfield Q into the 50(10) representations. It is useful to consider the 
decomposition of the 10 and 16 of 50(10) with respect to the subgroup SU{5) U{l)x, 
namely 

16 = 10i©5_3©l5, 10 = 5_2e52, (2.1) 

while the decomposition of the 16 follows from that of the 16. In order to avoid negative 
soft terms contributions (cf. Eq. (1.2)), all the MSSM matter superfields {q,u'^ ,e'^ ^d!^ ,€) 
must have same sign under U{l)x- This condition is fulfilled if 

g © u'^ © = lOi C 16 , (i^^ © ^ = 52 C 10 . (2.2) 

The MSSM Higgses are embedded in linear combinations of 10, 16 and 16. According to 
this discussion the TGM contribution to the bilinear soft masses is 



(rfif 



Q J tree 



2m?o Q = d^,i 



10 

(2.3) 



"^10 Q = e' 



-2mfo < (m^Jtrcc < 3mfo 
-3?nfo < i^'^hj tree < 2mfo 



where m^Q is a universal mass parameter. 

Gaugino masses are generated at the one-loop level as in standard gauge mediation. 
We call the common gaugino mass at Mq^ ■ 

The mechanism responsible for gaugino masses unavoidably generates a two-loop level 
contribution to the sfermion masses. The final expression at the GUT scale is given by 

= (^Q)trce + 2 7? CqM^^ , (2.4) 



^We are not taking into account here the efTect of having a hierarchical messenger spectrum, an analysis 
of this scenario will be presented in Ref. [6]. 



-3- 



where Cq is the total SM quadratic Casimir 



Q 


q u'^ dp £ hu ha 


Cq 


21/10 


8/5 7/5 9/10 3/5 9/10 9/10 



(2.5) 



The parameter rj > gives the relations between the two-loop contribution to the sfermion 
masses and that to the gaugino masses squared'^. The precise value of r] depends on the 
details of the messenger sector. For instance in standard (loop) gauge mediation with one 
messenger chiral superfield this parameter is precisely rj = 1/n, where n is the Dynkin 
index of the vector- like pair of messengers. In most of our analysis we will set r] = 1. 

In order to clearly discriminate TGM as the mechanism responsible for the soft scalar 
masses we will focus the attention on the regions of the parameter space where the values 
of mio and such that TGM is responsible for the leading contribution to the 

sfermion masses. We define the dominance of TGM by requiring that TGM contributes to 
the low energy value of each sfermion mass by an amount of at least 50%. Including also 
the running effects (cf. Appendix A for further details) this translates into the condition 

mio > (5.2 + 4.2r?)i/2^fi/2, (2.6) 

which, for ?7 = 1, reduces to mio ^3.1 Mi/2- 

On top of TTT-io and Mi^2 the other MSSM parameters relevant at low-energy are tan/3, 
fi, Bfi and the A-terms. Relating the //-term to SUSY breaking is a model-dependent 
issue'^. Here we will just fix fi and Bfi in such a way that they satisfy the electroweak 
symmetry breaking conditions. In addition we assume fi > 0. The ^d-terms are set to zero 
at the GUT scale. In general, since they do not arise at the tree level, they are expected 
to be smaller than the bilinear soft masses. 

In the case in which SUSY is broken only by the F-term responsible for sfermion 
masses, the gravitino mass is directly related to rhio by the relation [1] 

-3/2-15Gev(^) . (2.7) 

We stress again, due its cosmological relevance, that the magnitude of the gravitino mass 
is a peculiar prediction of TGM and its embedding into a GUT. 

Let us close this section with a couple of comments regarding the spectrum. The 
first one is about the nature of the NLSP, being a cosmologically relevant issue. It turns 
out that in most of the parameter space the NLSP is a Bino-like neutralino. This is 
easily understood because, according to Eq. (2.6), scalars are expected to be heavier than 
gauginos. Moreover, the running of gaugino masses (cf. Eq. (A.l)) is such that Mi < M2. 
The second comment regards the viability of the model in the light of the recent LHC 
exclusion limits [8]. The condition in Eq. (2.6) can be translated (cf. Eqs. (A.20)-(A.21) 

^In the notation of [7] this parameter can be identified with the ratio 77 — A|/Aq. 

*We just mention that TGM offers a new solution for the /i-problem, where the fi is also responsible for 
triggering the SUSY breaking [2]. 
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in Appendix A) in terms of the values of the physical stop mass rht and the gluino one rrig, 
yielding 

fht > 1.2 TeV ( \ , (2.8) 

V700 GeV/ ^ ^ 

which shows that a sizable part of the parameters space is still viable and testable at the 

LHC. 



3 Gravitino Dark Matter with R-parity 

In this section we investigate the cosmological aspects of TGM in the R-parity conserving 
case. The cosmology is deeply influenced by the behavior of the NLSP. Indeed, due to the 
presence of R-parity, it affects the DM relic density by decaying into gravitinos. Being the 
rate of such a decay Planck-suppressed, gravitinos will be produced after the freeze-out of 
the NLSP and potentially also after the onset of the BBN. We then refer to these gravitinos 
as non-thermal and their abundance is given by 

^3/2^ = — ^NLSph^ ~ 3 X 10 ^— — n^LSph^ , (3.1) 

' rriNLSP Ml/2 

where ^^nlsp is the expected NLSP relic density as if it were stable. 

ri^^ has to be added to the abundance J^^^/ of gravitinos produced by thermal pro- 
cesses in the early Universe, in order to match the WMAP-7 value Q.^uh'^ = 0.1123 it 
0.0035 [9]. The thermal component can be computed once the rates of the relevant pro- 
cesses are known as a function of the MSSM parameters (see e.g. [10]) and it depends 
on the reheating temperature T^h, i-e. the temperature which sets the beginning of the 
radiation domination era. 

The non-thermal component for the model in consideration is obtained by means of 
the scaling formula in Eq. (3.1) with i^NLSP computed through the numerical code DARK- 
SUSY [11]. The results are reported in the left panel of Fig. 2. The outcome is that the 
scenarios in which TGM is the dominant mechanism for the generation of sfermion masses 
(cf. the region in Fig. 2 above the blue line) are characterized by an overabundance of 
non-thermal gravitinos. 

This behaviour can be easily understood from the fact that Binos annihilate mainly 
into fermion pairs with a p-wave suppressed cross-section oc m'^o/mQ, where m^o is the 
lightest neutralino mass. In such a case the NLSP relic density can be estimated by 

W/.= »0.02xl0'(^)"(^)'. (3.2) 

Combining this result with Eq. (3.1) it is evident that the cosmological value of the DM 
relic density is largely overcome for mio > Mi/2i ^ predicted by TGM. 

The decay of the NLSP is also responsible for another important cosmological issue. 
Indeed the gravitino injection is accompanied with the production of SM particles which 
can trigger either electromagnetic or hadronic showers. In turn these may upset the BBN 
calculations of the light-element abundances. As shown in the right panel of Fig. 2 the 



- 5 - 




500 1000 1500 2000 500 1000 1500 2000 

M|,2 [GeV] M|/2 [GeV] 



Figure 2: Contours of the non-thermal component of the gravitino rehc density (left 
panel) and of the neutralino lifetime (right panel) in the plane {rhiQ, M1/2), with the 
other MSSM parameters fixed as explicitly said in the text. The black line represents the 
non-thermal relic density fit of the WMAP-7 value. The blue line represents the relation 
mio = 3.1 which sets the TGM dominance. 



lifetime of the neutralino is much bigger than the time of the onset of the BBN also for 
moderate values of the ratio rhio/Mi/2. 

We can easily understand this behaviour from the functional dependence of the decay 
rates. A Bino NLSP mainly decays either into a Z boson and a gravitino or into a photon 
and a gravitino and the relevant rates can be expressed as [12] 

nA^ZG)='^^, r(x?^7G) = ^^. (3.3) 
^ 487rM|, m2^2 ^SirMj, ra\i^ ^ ' 

In particular the NLSP can induce hadronic showers from the decay of the Z boson. These 
have the deepest impact on the BBN and hence suffer from the most severe bounds. The 
constraints are both on the lifetime and the abundance of the NLSP, and they basically 
exclude all neutralinos with a lifetime greater than 10~^ s [13]. We conclude that is not 
possible to obtain a viable cosmology in the TGM setup with R-parity conservation. 

It is also evident that the BBN bounds, together with the issue of the overabundance 
of gravitinos, can be evaded in the presence of some mechanism which suppresses the NLSP 
abundance. In this paper we will consider the case of a small amount of R-parity violation^. 

For completeness we mention that, by relaxing the condition Eq. (2.6), it is possible 
to realize scenarios with a viable cosmology. Here we report two examples. The first one 



^An alternative scenario could be the dilution of the NLSP abundance due to the entropy released by 
the decay of some heavy state before the onset of BBN [14, 15]. 
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is given by the enhancement of the s-channel annihilation of neutrahnos into a bottom 
pair mediated by the CP-odd Higgs, typicahy occuring at high values of tan/3. This 
enhancement is particularly strong in the region uia ~ 2m,-;^, where the annihilation cross 
section can also become resonant. As shown in Fig. 3, this scenario is realized only if one 
deviates sensibly from the relation Eq. (2.6) and hence in a region where standard (loop) 
gauge mediation is the dominant mechanism for the sfermion mass generation. 
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Figure 3: Left panel: WMAP-7 value of the non thermal contribution to the DM relic 
density, within a 1-a deviation. Blue dot-dashed lines represent the contours of the 
NLSP lifetime. Right panel: points of the left panel plot which give a non-thermal grav- 
itino relic density lower (red points) or within 1-a (green points) with respect to the 
WMAP-7 expectations. The blue and the violet lines represent, respectively, the BBN 
bound for a 100 GeV and 1 TeV decaying particle [13, 16]. 



The second example is the case in which the NLSP is a stau. This is possible for high 
values of tan/3 and r/ < 1. In this case the negative Yukawa corrections to the third family 
sfermions, originating at high tan/3, can drive the mass of the lightest stau below that of 
the lightest neutralino. The stau mainly decays into a gravitino and a tau lepton with the 
following rate [12] 

r(f^TG)=—\-^^, (3.4) 

without producing hadronic showers. In this case the strongest BBN bounds are given by 
the formation of bound states with Helium nuclei which can alter the primordial abundance 
of Lithium. This process is referred as catalyzed BBN [17, 18] and implies an upper bound 
of around 5 x 10'^ s on the stau lifetime^. Given the rate Eq. (3.4) this translates into the 
requirement of a stau mass greater than 200 300 GeV. 

®It should be mentioned that in settings with substantial left-right mixing of the stau mass eigenstates 
the BBN bounds can be evaded even for higher lifetimes [19]. 
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4 Gravitino Dark Matter without R-parity 

The outcome of Sect. 3 is that the configurations of TGM dominance, rhiQ > 3.1Mi/2i 
cannot fulfill the cosmological bounds regarding the DM relic density and the BBN. On 
the contrary these bounds can be evaded in presence of a small amount of R-parity breaking 
(see e.g. Ref. [20]). Indeed in such a case the NLSP can decay only into SM particles before 
the onset of the BBN, without contributing to the DM relic density in the form of non- 
thermal gravitinos. At the same time the thermally produced gravitinos should be stable 
enough in order to reproduce the correct DM relic density. Then, since the NLSP lifetime 
does not depend anymore on the gravitino mass, a hierarchy between sfermion and gaugino 
masses, as expected in TGM, is in principle achievable. 

In the following we will present and analyze an 50(10) TGM model with explicit R- 
parity violation. In particular we will show that applying the bounds from proton decay the 
GUT structure of the theory guarantees the suppression of all the trilinear RPV couplings. 
This motivates the restriction to only bilinear RPV operators in the cosmological analysis. 

4.1 An R-parity violating 50 (10) model 

The R-parity is a Z2 symmetry which distinguishes the SM fields from their super-partners. 
By assigning a positive charge to the former and negative one to the latter it provides for 
instance the stability of the LSP. It is useful to rephrase the R-parity in a slightly different 
language [21] 

Rp = {-f(^-^)+^\ (4.1) 

where s is the spin quantum number. As long as the angular momentum is conserved, the 
R-parity is essentially equivalent to a Z2 Matter (M)-parity, defined as 

Mp = (-)3(s-^) . (4.2) 

Eq. (4.2) suggests that in theories in which B — L \s gauged the M-parity can be viewed as 
a discrete subgroup of a local U (1) _b_l, thus providing a possible link between the amount 
of R-parity violation and the B — L breaking scale. However, though one of the S'0( 10) 
generators can be identified with B — L when acting on the spinorial representation 16 
this does not apply to TGM because of the peculiar embedding of the SM fermions into 
16^0 IOf (cf. Eq. (2.2)). 

Thus the 50(10) gauge symmetry does not protect against the appearance of R-parity 
violating interactions and the simplest way to forbid them is to impose an M-parity external 
to 50(10) which distinguishes the matter superfields (with negative charge) from the Higgs 
ones (carrying positive charge) [1]. 

Our approach will be that of relaxing the presence of this extra M-parity and, by 
considering all the operators compatible with the 50(10) symmetry, we will prove the 
existence of R-parity violation in the low-energy effective theory. In addition we will 
assume that the R-parity violating operators are suppressed when compared to their R- 
parity conserving counterparts, in order to avoid an unacceptable amount of lepton and 
baryon number violation. Addressing the issue of the origin of such a small amount of 
R-parity breaking is anyway beyond the scopes of this work. 
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For definiteness let us focus on a supersymmetric 50(10) model featuring the follow- 
ing minimal set of Higgs representations: 54// © A5h © ISjy © IQh © lOj/^. As shown 
in Appendix B this field content is sufficient in order to break 50(10) down to the SM at 
the renormalizable level (cf. Appendix B.l) and to give mass to the SM fermions (cf. Ap- 
pendix B.2). As already mentioned the MSSM matter superfields span over three copies of 
16i? © IOf in such a way that they are embedded in the SU{5) representations 10 D 16/? 
and 5 D Wp (cf- Eq. (2.2)). The conditions to be fulfilled in order to obtain such a "pure" 
embedding are detailed in Appendix B.2. 

The superpotential can be schematically written as 

W = Wh + Wy + 6Wrpv , (4.3) 
where Wh and Wy are the Higgs and the Yukawa components 

Wh = (^54 + f?5445// + A5454h) 54|/ + /U45 45|/ + (^no + Aio54//) lOjj 

+ (/ii6 + Ai645i/) 16hT6h + Ai6-iol6^10// + Ai6-iol6^10// , (4.4) 
Wy = Yi^lG^lG^Wn + yHigUo'fIQh + + ??'^45h + X'^BAh) 10^10^, , (4.5) 
while 6Wrpv is the R-parity violating piece 

6Wrpv = (/ilo + ^io45H + Aio54//) 10^10// + + ~X\q45h) 16>T6i/ 

+ p'IG^WhWh + a^lOVl6//16// +^'10VT6//T6// + A^^'=165.16],10^ . (4.6) 

Notice that without M-parity the separation between the F and the H superfields is some- 
how artificial. However, since we consider 5Wppv as a perturbation, we can still retain Wh 
responsible for the symmetry breaking and Wy for the (charged fermions) Yukawa sector. 
On the other hand the situation about neutrino masses is subtler, being RPV potentially 
responsible for sizable contributions to them. We will comment later on the generation of 
neutrino masses in our model. 

In Appendix B.3 we provide an existence proof of the R-parity violating operators in the 
MSSM effective theory. In order to obtain the low-energy superpotential one has to project 
the operators of Eq. (4.6) on the representations containing the MSSM fields. Here we just 
report the results of this operation leaving most of the technical details in Appendix B.3. 

Bilinear i?-parity violation in the effective superpotential is induceded by operators 
containing just one F superfield in Eq. (4.6), leading to 

w'J/yD ^l^eihu, (4.7) 

where the expression of Hi in terms of the original couplings is given in Eq. (B.25). 

Notice that in the "pure" matter embedding of TGM only some of the operators have 
projections on the MSSM fields. In particular the trilinear operator relative to the coupling 

^To be complete one should also add a 16'i] © IG'h representation which is responsible for the SUSY 
breaking [1]. For simplicity we will carry on our analysis in the supersymmetric limit, assuming that the 
conclusions regarding the gauge symmetry breaking and the fermion mass spectrum are only marginally 
affected by the SUSY breaking sector. 
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Jyijk does not contribute to the effective tlieory. On the other hand the phenomenological 
viabiiity of the model, within the minimal choice of representations at hand, requires the 
presence of non-renormalizable operators (cf. again Appendix B). By relaxing renormaliz- 
ability there is an additional source of R-parity violation given by the operator 

10^^10].16^ (16^> D Ajf Id A + A^f,^ dii.q, + X^Jl^ d^d^vfi , (4.8) 



Mp 



where 



\jk — l^^jk — \jk — ^[JZ = ■ ^ 



Notice that the relation in Eq. (4.9) gives a correlation between the baryon (A") and lepton 
number (A, A') violating couplings. 

We should also mention that A" receives an additional contribution when combining 
the bilinear operators in Eq. (4.6) with the Yukawa ones. This is obtained by projecting 
the Higgs fields on the heavy triplet components and integrating them out. This last 
contribution, labeled A^, is shown explicitly in Eq. (B.29). 

In the end the structure of the induced superpotential in the MSSM effective theory 
is given by: 

= /i. ^^K + idjel + a:,., dlljQk + A'/., dtdy, , (4.10) 



where 



A = A, A' = 2A, A" = A + A^. (4.11) 



The strongest constraints on the R-parity violating interactions are due to proton decay. In 
particular this enforces severe bounds on the products of couplings AA", A' A" (cf. e.g. [22] 
for an exhaustive list). For kinematical reasons the most stringent ones apply to products 
involving dominantly the first two light generation indices. For instance the process p — >• 
TT^e'^ severely constrains the product 

A;nA?„<10-(^)\ (4.12) 

with k = 2,3 and rh being the sfermion mass scale. Though the extension of the analysis 
at the one-loop level sets weaker bounds for the couplings relative to the second and third 
generation [23], in our setup the structure of the trilinears is constrained by the GUT 
symmetry which implies much stronger bounds with respect to the general case. The most 
conservative bound on all the R-parity violating trilinear couplings in the presence of a 
GUT relation such as that in Eq. (4.9) is given by [24] 

Barring extremely accurate cancellations between the two unrelated components A and Xt 
in the expression for A" (cf. Eq. (4.11)), the bound in Eq. (4.13) is automatically translated 
onto A and A'. 
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The bounds on the trilinear RPV couphngs just derived are very strong, making them 
harmless for the cosmological analysis. In light of this result, a motivated setup for the cos- 
mological analysis is bilinear R-parity violation. Then the R-parity violating superpotential 
simply reads 

w''JJiy = ^l^^iK. (4.14) 

In the effective theory one also expects R-parity violating couplings in the soft scalar 
potential, depending on the details of the SUSY breaking sector. As we will show in the 
next section, the main cosmological constraints apply to the bilinear soft terms 

yRPV = Bi hhu + rhl^i^ h\ii + h.c. . (4.15) 
4.2 Cosmological analysis 

Without R-parity, and thus lepton number conservation, there is no a priori distinction 
between and li. This is reflected by the fact that Eqs. (4.14)-(4.15) induce a non- 
vanishing VEV for the sneutrino fields once electroweak symmetry is broken [22, 25]. 

Since and li have the same quantum numbers one can always operate a redefinition 
of these fields through a unitary transformation. Typically one can use this transformation 
in order to rotate away either the sneutrino VEVs Vi or the bilinear couplings /ij. 

Since we consider Eq. (4.14) as a perturbation of the R-parity conserving theory, 
implying in particular /Xj <C /x, it is convenient to define the following linear transformation 
on the superfields 

hd^ hd = hd + eiii , ii ^ ii = ii - eihd , (4.16) 

with = fJ-i/fJ., which rotates away the bilinear term in Eq. (4.14) up to 0{ef) corrections. 
The expressions of the R-parity violating couplings in the hatted basis can be found for 
instance in Ref. [22]. We just point out that the transformation in Eq. (4.16) induces 
trilinear lepton violating couplings Xijk, X'^j^. of the form 

\ijk = -{Ye),k^j + {Ye)^f,e,, Kjk = -{yd)ik^j, (4.17) 

where Ye and Yd represent, respectively, the SM Yukawas of the charged-leptons and the 
down-quarks. In the hatted basis the VEVs of the sneutrinos are given by [20] 

Bi tan B + mi n 

V. = -Uhd) = -^T-T^T^^(M . (4.18) 
m|. -I- iM|cos2/3 

Then we can express in terms of the original parameters obtaining, at the leading order 
in ei, 

(5i - ei-B) tan /3 m|^^^ (r 



. , . (4.19) 

m^, + cos 2/3 ^ ' 

Given the model dependence of the soft terms Bi and rhh^i^ , RPV is described at low-energy 
by the six parameters and ej. On the other hand, by inspecting Eq. (4.19) and Eq. (2.3) 
it turns out that, barring cancellations, the parameters are at least of the order of ej. 

In what follows we will analyze the bilinear RPV model introduced in the previous 
section and identify the range of viability of the R-parity violating couplings. The analysis 
is organized according to the following points: 
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• NLSP lifetime and BBN. 



• Gravitino lifetime and cosmic rays. 

• Neutrino masses through RPV. 

• Gravitino relic density and thermal leptogenesis. 
NLSP lifetime and BBN 

The main motivation for the introduction of R-parity violation is to restore the agreement 
between the decay of the NLSP and BBN. In our setup the most efficient processes are 
induced by the sneutrino VEVs which mixes the Z (VF) boson with a neutrino (charged- 
lepton) and a neutralino [26, 27]. 

A Bino can decay into a W boson and a charged-lepton or into a Z boson and a 
neutrino with the typical rates [28] 



where ^ = vCi+W+W- '^^^ NLSP may also decay into three fermions by means of the 
couplings A and A' in Eq. (4.17) with a typical rate of the form 

Ta-body = ino/i • (4.22) 

l(J24 7r'^ ttIq 

The same expression, divided by a factor of three, holds for the rates involving the coupling 
A. However the 3-body processes are highly suppressed with respect to the 2-body decay 

BRC, - body/2 - body) » 0^9 x 10- f 1 V f 7^ V7 ' . (4^23) 



\ J \l TeV J \lbO GeV 

In writing this expression we used Eq. (4.17) and assumed ej ~ e for i = 1, 2, 3. Hence the 
3-body processes can be neglected for a typical TGM spectrum, barring cancellations in 
the sneutrino VEVs. The NLSP lifetime is then determined by Eqs. (4.20)-(4.21) 

™LSP^.-.oa, - 0,02 s (^^) ' (^) ' (^) " . (4,24) 



which satisfies the BBN bounds, as reproduced in the right panel of Fig. 3 for tnlsp ^ 10 ^ 
s, thus implying ^ > io~(^'^~^^) depending on m^o. This is actually a rather conservative 
bound; lower values of ^ are allowed depending on the NLSP abundance®. 



* On the other hand the given bound takes also into account the fact that a population of NLSP, poten- 
tially dangerous for BBN, could survive despite the low branching ratio of 3-body decays. From Eq. (4.23) 
we see that this scenario does not occur for low/moderate values of tan P (cf. also the lines in the second 
panel of Fig. 3). Eventually it could be necessary to assume an even more conservative limit ^ > 10~^. 
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Gravitino lifetime and cosmic rays 

The amount of R-parity violation is also constrained from above since the DM is not 
stable anymore. As it will be evident from the expressions below the gravitino is stable 
on cosmological time-scales, being its decay rate doubly suppressed both by the R-parity 
violating couplings and the Planck mass. On the other hand the small portion of the 
decaying gravitinos is able to leave an imprint on the cosmic ray spectrum. 

In our setup the main decay channel of the gravitino is into a neutrino and a photon [29] 

Then the associated lifetime can be estimated by^ 

, 7.3 X 10- s f^^V (^) " f « ) " . (4.28) 



10 J V300 GeVy V15GeV/ VlO 

This process is expected to leave an imprint on the cosmic gamma ray spectrum in the 
form of an approximately monochromatic line at an energy depending on the gravitino 
mass. For higher values of M1/2 and 7713/2 the 3-body processes mediated by off-shell gauge 
bosons become also important and eventually dominant [31, 32], implying the presence of 
an additional continuos component in the gamma ray spectrum. This kind of signals have 
been the subject of dedicated searches performed by the Fermi Large Area Telescope (LAT) 
collaboration [33, 34]. Since none of the expected excesses in the gamma ray spectrum have 
been detected so far one obtains a lower bound on the gravitino lifetime. 

According to Refs. [31, 35] the lower bound on the gravitino lifetime is approximatively 
of 10^^~^^ s for gravitino masses in the range 10-^80 GeV and ™ the range 100-^1000 
GeV. This translates into an upper bound for ^ of about ^ < 10"'-^"'^^ 

For definiteness we mention that for the central values of M1/2 and 7713/2 in Eq. (4.28), 
which roughly correspond to the LHC bound quoted in Sect. 2, the limit on the gravitino 
lifetime is 10^^ s which is satisfied for C ^ 3 x 10~^. 

Neutrino masses trough RPV 

Summing up the results obtained until now, the outcome of the cosmological analysis is 
that the RPV coupling ^ must lie in the window 10"^^ < ^ < 10^^. R-parity violation 
can have, however, a wider impact on the phenomenology being a potential source of 



^The trilinear couplings A and A' yield a negligible contribution to the gravitino lifetime. Indeed the 
typical rates are [30] 

A^ ^3/2 - -I - 

r3/2,3-bod.= ^g^32^3 , A = 3A,A (4.26) 

which lead to the following lifetime 
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neutrino masses besides its cosmological role^*^. Indeed the VEVs of the sneutrinos induce 
a mixing between neutrinos and neutrahnos leading to a contribution to neutrino masses 
after integrating out the neutrahnos [25]. The result is a rank-one matrix hence giving 
mass to just one neutrino. Explicitly one finds [37] 

m.3 = M|eW/3( ff' , - ^ sin 2/3)"' . (4.29) 

The complete neutrino spectrum can be then reconstructed by including also the one-loop 
contributions which depend on A, A' and the soft parameters Bi (see for instance [38, 39] 
for a detailed computation of neutrino masses at the one-loop level). Barring cancellations 
in ^ the neutrino spectrum turns out to be hierarchical, hence we can impose the relation 
m^.j ~ a/ Am^^-jjj, thus implying 

Additional constraints on the parameters and can be imposed once other observables, 
such as Arrig^j and the mixing angles, are taken into account. As it is evident from Eq. (4.30) 
the outcome of the analysis of cosmic ray bounds (^ < 10~^) implies a suppression of 
neutrino masses induced by RPV well below the experimental constraints. This result 
is in agreement with the recent analysis of Ref. [40] which rules out bilinear RPV as 
the mechanism responsible for neutrino masses for gravitinos heavier than 1 GeV. Such 
low vales of the gravitino mass are not achievable if TGM is the dominant mechanism 
originating sfermion masses (cf. Eq. (2.7)). In our model, neutrino masses have to be 
generated by means of another mechanism (see Appendix B.2 for the discussion on neutrino 
masses) . 

Gravitino relic density and thermal leptogenesis 

At this point we turn to the relic density of DM. For the range 10^(^'^~^i) ^ C ^ 10^ 
the NLSP decays only into SM particles at a much faster rate with respect to the decay 
into gravitinos. The branching ratio between the R-parity conserving decays (into SM 
particles) and the R-parity violating ones (into SM particles and gravitinos), 

indicates that the DM relic density is completely determined by its thermal component. 
This can be computed as a function of the reheating temperature Tjih in thermal field 
theory, according to the following analytic expression [41] 

i'l'DM'^ — 



10 GeV/ V 109 GeV 



.)E.;.?(«a.«(i.^)i..(^). (4.32) 



^"in this respect it has been pointed out that the MSSM without R-parity could be a complete theory 
of the low-energy phenomena including neutrino masses, and it could even fit the cosmic ray anomalies of 
PAMELA and Fermi-LAT in terms of (decaying) gravitino DM [36] . 
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where r = 1, 2, 3 and the sum runs over the three components U{1)y, SU{2)i and SU (3)c 
of the SM gauge group. The values of the coefficients k^-, y'^ and br can be found for instance 
in Ref. [42]. 

As shown in Fig. 4 the DM relic density matches the cosmological value, depending 
on the values of and mio, for reheating temperatures up to ~ 10^ GeV. By estimat- 
ing Eq. (4.32) in the following way 

and considering the relation in Eq. (2.7), we can see that the increase in the contribution 
of TGM to sfermion masses (with respect to standard gauge mediation) coincides with an 
increase of the reheating temperature needed to fit the cosmological value of the DM relic 
density. Moreover values of the reheating temperature of the order of 10^ GeV are welcome 
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Figure 4: Contours of the gravitino relic density in the plane (Mi/25"^io) computed, 
according to Eq. (4.32), for the values of the reheating temperature reported. The green 
dashed vertical line represents the value of 700 GeV for the gluino mass while the black 
line represents the TGM dominance relation in Eq. (2.6). 

in the standard thermal leptogenesis scenarios [43]. We stress that this result is obtained 
thanks to the natural prediction of TGM of a gravitino mass of the order of 10 GeV. On the 
contrary theories based on standard (loop) gauge mediation predicts much lower reheating 
temperatures due to lower gravitino masses. By rescaling Eq. (4.33), in terms of mio and 
the gluino mass 




we see that a reheating temperature of the order of 10^ GeV requires sfermion masses in 
the multi TeV range. 



We conclude by mentioning that, in presence of a leptogenesis mechanism, the theory 
is subject to additional constraints on the amount of RPV. Indeed baryon and lepton 
number violating interactions due to RPV could erase the B — L asymmetry generated 
by leptogenesis. This can be avoided by requiring that the dangerous processes are not 
efficient, i.e. Trpv < when the asymmetry is generated. The general expression for these 
rates has been studied in [44, 45]. In case of bilinear R-parity violation the baryogenesis 
bound can be summarized by the condition [20] 



e < 10^ 



tan /3 
10 



(4.35) 



which implies a similar bound on ^, barring cancellations in Eq. (4.19). 



The bounds emerging from the cosmological analysis are collectively summarized in Ta- 
ble 1. 



Observable 


Bound 


References 


Proton decay 


A, A', A" < 10-10 


[24] 


BBN (NLSP hfetime) 


e > lo-(io-ii) 


[46] 


Cosmic rays 


i < 10-(6-8) 


[31, 35] 


Neutrino masses 


^ ~ 10-5 


[22] 


Baryogenesis 


e < 10-6 


[44, 45] 



Table 1: Summary of the bounds on the R-parity violating couplings. 



5 Conclusions 

In this paper we have discussed the impact of cosmology on the 50(10) TGM model. The 
new gauge mediation mechanism introduced in [1] guarantees the flavor universality of the 
sfermion masses and peculiar predictions for the supersymmetric spectrum at the LHC. 
The messenger scale is fixed at the GUT scale and this gives also a prediction for the 
value of the gravitino mass to be of 0(10 GeV) by assuming the sfermion masses in the 
TeV range. Another consequence of the model, under the assumption of TGM dominance 
(cf. Eq. (2.6)), is the fact that the NLSP is a Bino-like neutralino. 

This scenario is highly disfavored if R-parity is conserved, being the DM overproduced 
by the NLSP decays. Moreover the gravitino production is accompanied by hadronic and/or 
electromagnetic showers which spoil the predictions of BBN. 

On the contrary the TGM model is naturally feasible in presence of a small amount of 
R-parity violation. Furthermore the GUT structure of the theory and the constraints from 
the proton decay motivates the restriction to a scenario of bilinear R-parity violation which 
allows to describe the relevant phenomenology in terms of a limited number of parameters. 
Given the correct amount of RPV, the NLSP is allowed to decay much before the onset of 
BBN without overproducing gravitinos. The latter, at the same time, remain stable over 
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cosmological times, being their decay rate doubly suppressed both by the R-parity violating 
couplings and by the Planck scale. Interestingly the small amount of decays which can take 
place at the present time is in principle detectable by the current cosmic ray experiments 
such as FERMI-LAT. 

The same RPV couplings responsible for the NSLP and gravitino decay could be at the 
origin of neutrino masses. However the value required for these couplings is not compatible 
with the bounds from cosmic rays, at least for the gravitino masses predicted by TGM. 

On the other hand a gravitino with mass of 0(10 GeV), combined with a natural 
TGM spectrum with sfermions in the multi TeV region, allows the DM relic density to 
match the cosmological value for a reheating temperature of 0(10^ GeV) relevant for 
leptogenesis. This improves the situation with respect to the case of the standard (loop) 
gauge mediation, where the gravitino mass can be at most of 0(1 GeV) implying a lower 
reheating temperature. 
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A TGM dominance for sfermion masses 

In this appendix we perform an analytical study of the one-loop renormalization group 
equations (RGBs) for the relevant MSSM parameters in order to investigate the nature 
of the NLSP and to provide a criterium for the determination of the TGM dominance to 
sfermion masses. Defining t = log(/x//xo), the running of gaugino masses for a = 1,2,3 
reads 

Ma(i)=Mi/2^, (A.l) 

% 

where go and M1/2 are evaluated at the GUT scale and ga{t) obeys the RGE 

9a\t)=g^'-^t, (A.2) 

with ba = (33/5, 1, —3) being the one-loop MSSM beta-functions. Assuming flavor univer- 
sality at the GUT scale and neglecting the contributions from the Yukawa couplings and 
the A-terms^^, the one-loop RGEs for the scalar soft masses can be written as [47] 

167r2^m| = - J2 8Ca{Q)glM! + ^YQglS, (A.3) 



^^We can safely assume that the Yukawa terms give no sizable contribution to the sfermions of the first 
two families. This is not necessarily the case for the third family, but since we are interested in a simple 
analytical estimate of the running effects we stick anyway to this approximation. 
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where Ca{Q) is the quadratic Casimir relative to the representation Q and to the gauge 
group a, and the factor S reads 

S = '^L- "^L + 3 (mg - - 2ml, + ^% + rh%) ■ (A.4) 
Then, by combining Eq. (A. 3) and Eq. (A.4) we obtain 

l67r'p=fglS, (A.5) 

whose integration yields 

S{t) = , (A.6) 

% 

where Sq is evaluated at the GUT scale. Given the SO{W) embedding of the matter 
superfields in Eq. (2.3) one gets Sq = 5(0) = m^^(O) — ?7i^^(0). 

We can now integrate Eq. (A. 3) with the help of Eq. (A.l), Eq. (A. 2), Eq. (A.6) and 
the GUT boundary conditions in Eq. (2.4), obtaining 

mlit) = Aq mjo + {BQ{t) + 2 Cq r?) Af,^ + Dgit) So , (A.7) 

where Aq is equal to 1 for z = q, u^, e'^ and 2 for i = ^, (F (cf. Eq. (2.3)) and Cq is the total 
SM quadratic Casimir relative to the representation Q (cf. Eq. (2.5)). The coefficients Bq 
and Dq parametrize the effects of the running and can be analytically expressed as 

BQ{t)= Y: 2C.(g)l(l-^) and i^Q(t) = ^ - l) • (A.8) 

a=l,2,3 °" ^ ^0 ^ V 5o / 

Evaluating Eq. (A.l) and Eq. (A.7) at the TeV scale one obtains 



Ml 


-0.4 Ml/2, 


(A.9) 




^0.8 Ml/2, 


(A.IO) 


Ms 


-2.5 Ml/2, 


(A.ll) 


ml 


~ + (5.2 + 4.2 r/) M^^^ - 0.009 Sq , 


(A.12) 


m\c 


~ m?o + (4-8 + 3.2 r/) M^^^ + 0.03 So , 


(A.13) 


m% 


~ m?o + (0-1 + 3-2 ?7) m2/2 - 0.05^0 , 


(A.14) 


m% 


~ 2 m?o + (4.7 + 2.8 r/) M^/^ - 0.02 Sq , 


(A.15) 


mj 


~ 2 mfo + (0.5 + 1.2 ??) M2/2 + 0.03 5o . 


(A.16) 



The presence of the contributions due to M1/2 makes TGM the leading mechanism for 
sfermion masses only in some portions of the MSSM parameter space. We can define oper- 
atively the dominance of TGM by requiring that at least the 50% of the low-energy value 
of the sfermion masses is due to TGM. Neglecting 5o the worst case scenario in Eq. (A.12) 
translates into 

m2o>(5.2 + 4.2 77)M2/2, (A.17) 
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that is 

mio> 3.1 Ml/2, (A.18) 

for 77 = 1. From this relation it is evident that the NLSP is always the lightest gaugino if 
TGM is the dominant mechanism generating sfermion masses. 

Given M1/2 — 0.4 M3 from Eq. (A. 11) and taking M3 w irig, we arrive to the relation 

mio>1.2mg, (A.19) 

that, if substituted into Eqs. (A.12)-(A.16), yields the following bounds on the sfermion 
masses as functions of the gluino mass nig-. 



B Details of the 50(10) model 

In this Appendix we give the details of the 50(10) model presented in Sect. 4.1. For later 
convenience let us set the following notation for the SM components of the 50(10) fields 
relevant for the Yukawa sector 

16f = {D" e L)5_^ e (n^ e g e e^)io+, © (z^^)i+5 (B.l) 
IOf = {D ® L%_, ® {d'' (B i)^^^ (B.2) 
16h = (Tf © hf)-,_,^ © (. . .)io+, © (. . .)i+, (B.3) 

TEh = (tJ^© /i?)5+3 ® (• • Oto^, © (• • Oi-s (B.4) 
Wh = {Tl' © 05_, © (ri° © hf)-,^^ (B.5) 

where a self-explanatory SM notation is employed and the outer subscripts label the 
SU{5) (g) U{l)x origin. The SU(2)l doublets decompose as g = {u ® d), i = (z^ © e), 
L = {N(B E), = {E'^ © iV^), K = {K © hi) and hd = {h^^ © h^). 

B.l Symmetry breaking and doublet-triplet splitting 

The set of Higgs fields 54// © 45// © 16// © 16// is sufficient in order to achieve a renormal- 
izable^^ breaking of SO(IO) down to the SM (see e.g. Ref. [52] for the study of the vacuum 



^^With only 45h © I&h ® 16i/ at play the requirement of a supersymmetric vacuum at the GUT scale 
is such that the little group is 5[/(5) [48-50]. In order to reach the SM gauge group one can either relax 
renormalizability [49] or add a 54£f [48, 50]. Since the first option introduces a delicate interplay between 
the GUT and the Planck scale which may be an issue for unification and proton decay (see e.g. Ref. [51]), 
we choose the second option. 
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patterns). In particular, the SM gauge group is obtained as the intersection of the httle 
groups preserved by the following VEVs: 



(5%) = sU{A)c ® SU{2)l ® SU{2)r , (B.6) 

(45i^)B_^ = y^-L SU{3)c SU{2)l SU{2)n (E> U{1)b-l , (B.7) 

(45,^)^ = V^' SU{A)c SU{2)l ® ?7(1)r , (B.8) 

{16h) = V^^ SU{5). (B.9) 

With the minimal set of Higgs representations at hand we can explicitly check the feasibility 
of the doublet-triplet (DT) splitting. To this end we compute the mass matrices for the 
doublets (Md) and the triplets (Mr)- From Wh in Eq. (4.4) we get 



Mt={-^ „ - _ , (B.ll) 



defined, respectively, on the basis {hi^, h]^){hf , hf) and (T^^o, T„16)(T]°, T}^). The 
relevant Clebsch-Gordan coefficients can be found for instance in Ref. [53]. 

Two light Higgs doublets, hu and hd, can be obtained by imposing the minimal fine- 
tuning condition detMo ~ in Eq. (B.IO), while leaving at the same time the triplets at 
the GUT scale (cf. Eq. (B.ll)). Working for simplicity in the real approximation the light 
components read 

hu = cos Ouhl^ + sin 9uhl^ , = cos 9dhf + sin Oahf , (B.12) 

where, taking into account the minimal fine-tuning condition, 9^^^, are fixed in the following 
way in terms of the superpotential parameters 

tang^ = - .^45 ' tangrf = - . (B.13) 

Notice that in general 0„ ^ 0^. In particular, the projection of v}^'^^ and v^^'^^ on the 
electroweak VEVs Vu = {hu) and Vd = {ha) is 

= Vu cos 0u , vl^ = Vu sin Ou , vf = Vd cos 6d , v}^ = Vd sin 0d. (B . 14) 

Worth of a comment is the fact that the natural (without fine-tuning) implementation of 
the DT splitting requires the introduction of additional representations. A solution along 
these lines, in the context of an 50(10) model of TGM, has been put forward in Ref. [2]. 

B.2 Yukawa sector in the pure embedding 

Let us turn now to the Yukawa sector of the model. The fiavor structure of supersymmetric 
50(10) GUTs with extended matter sector {16f © IOf) has been extensively studied in 
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Refs. [54, 55]. On the other hand the mechanism of TGM requires a pecuUar embedding of 
the MSSM fields which must fit into SU{5) representations with positive X-charge where 
SU (5) (8) U {l)x C SO{10). This is needed in order to guarantee positive sfermion masses^^ 
(cf. e.g. Eq. (1.2)). 

Such an embedding is expUcitly shown in Eqs. (B.1)-(B.2), with the lower-case fields 
{q, u"^, d'^, £, e^, u^) labeling the MSSM degrees of freedom. In order for this to work one 
has to ensure that the vector-like pairs D'^ (B L and D ® (cf. again Eqs. (B.1)-(B.2)) 
pick up a super-heavy mass term, thus decoupling from the low-energy spectrum. 

After the electroweak symmetry breaking the fields with the same unbroken quan- 
tum numbers mix among themselves. As far as the charged fermions are concerned the 
superpotential Wy in Eq. (4.5) yields the following mass matrices 

defined respectively on the basis {u){u'^), {d, D){D^, d*^) and {E, e)(e'^, £'^). We also 
defined (see e.g. Refs. [54, 55]) 

Ma = Mio + r?yJ5_i - ^Ay5^ (B.16) 



^15 

/i 
2 V 5 



MA = Mio-r?yi^ + i^/fAy5\ (B.17) 



where Mio and A {rj) are symmetric (antisymmetric) matrices in flavor space^^. Thus, by 
inspecting Eq. (B.15), the decoupling of the vector-like pairs D^(BL and D®L^ is achieved 
by requiring 

Ma = Ma = 0. (B.18) 

We call this the pure embedding condition which gives the desired embedding up to 
Vd/y^^ ^ 1 corrections. Given the symmetry properties of the matrices Mio, rj and A 
and the need to keep the VEVs of the 45// and 54// switched on for the S'O(IO) symmetry 
breaking, the pure embedding condition translates into Mio = = A = 0. 

Let us consider now neutrino masses. In this case Eq. (4.5) is responsible for the 
following Major ana mass matrix 



M^ 



(QY^qv^^ YiqV^^\ 
■ YiQvf 
Xw+ Ma 
\. ■ . Xw^ J 



(B.19) 



deflned on the symmetric basis (A^, u'^ji', N^)- In Eq. (B.19) w± = ((1, 3, ±1)54^) denotes 
a pair of VEVs induced by Wh- The contribution to neutrino masses due to the VEV of 
the scalar triplets goes under the name of type-II seesaw. 



Strictly speaking what one has to require from a phenomenological point of view is that possible 
negative contribution to sfermion masses, originating from a non-pure embedding, are anyway subleading 
with respect to the positive ones [2]. For simplicity we stick here to the pure embedding limit. 
^^The reason being simply because 10 ® 10 = Is ® 4:5a © 54s. 
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In order to estimate the order of magnitude of the induced VEVs let us consider the 
fohowing piece of superpotential evaluated on the vacuum 



{Wh) D /i54 {{V^^? + w+w.) + A54 {w+w.V''^) + Alo {w+{vff + wMu?) (B.20) 
and require the F-term conditions Fyj^ = 0, which gives 

However, as explained above, the pure embedding condition requires A = so that the 
type-II seesaw contribution to light neutrino masses vanishes (cf. Eq. (B.19)). 

Sticking to a pure embedding one can still invoke non-renormalizable operators in order 
to give a mass to neutrinos through a standard type-I seesaw mechanism^^. Consider for 
instance the following Planck-suppressed operators 

and ^ 

IVlp iVlp 
They contribute to the light neutrino mass matrix after integrating out v'^, yielding 

mi = {YdY-^yS) (sin^ 0u?^ ~ {YdY-^yS) (sin^J^ _ ^^ ^4) 

This value is naturally too small and requires a fine-tuning in the Yukawa structure in 
order to restore the agreement with the experimental data. 

One should also keep in mind that in an R-parity breaking scenario there are new 
lepton number violating operators which contribute to neutrino masses as well. However 
the issue of neutrino masses with R-parity violation is tightly correlated with cosmology 
and it turns out that the size of the RPV couplings needed by neutrino masses leads, for 
the range of gravitino masses expected in TGM, to an unacceptable decay rate of gravitinos 
in view of the recent bounds on cosmic rays (cf. Sect. 4.2). 

The bottom line about neutrino masses is that they are naturally too small in the 
minimal SO{10) model in consideration, though it is always possible to fit them with a 
standard type-I seesaw mechanism (cf. Eq. (B.24)) due to the presence of unknown Yukawa 
structures which are not correlated with the charged fermion sector. 

On the other hand it is also easy to understand that by introducing additional repre- 
sentations in the game one can fit neutrino masses without too much fine-tuning. An inter- 
esting possibility is the introduction of a 54^ that couples to with strength 7. If the 
54^ does not develop a GUT scale VEV then the pure embedding condition in Eq. (B.18) 



^^Notice that in the pure embedding one has the 5(7(5) relation Md — MJ (cf. Eq. (B.15)) which is 
phenomenologically unacceptable. In this respect the presence of non-renormalizable operators is welcome 
in order to unlock that relation. 
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is automatically fulfilled with 7 7^ 0, yielding a type-II seesaw contribution to neutrino 
masses [56, 57]. The latter also provides an interesting leptogenesis mechanism based on 
the out-of-equilibrium decay of the Higgs triplets. On the other hand it has been pointed 
out in Ref. [56] that this mechanism requires an high reheating temperature of at least 
0(10^^ GeV), while the cosmological value of the DM relic density is fitted by reheat- 
ing temperatures pointing towards a standard thermal leptogenesis scenario based on the 
type-I seesaw (see Sect. 4.2 for more details). 

B.3 Origin of the R-parity violating operators 

This last section is devoted to the derivation of the R-parity violating operators in the 
effective MSSM theory. Starting from 6Wrpv in Eq. (4.6) and by projecting the 5*0(10) 
representations onto the light components (cf. Eqs. (B.1)-(B.2)) one finds: 

• A bilinear operator of the type Hiiihu, where 



= COS 9u (jl\o - vloV^'' + ^^f Alo^^') + sin^^^V^^ . (B.25) 

• Two bilinear operators of the type fi^^d'^T^^ and n^^^d'^T^^ , where 

^}p^ = j2^^ + f^l^Vit^-^^~X\oV'^ and # ' = ^V^^ (B.26) 

The triplet bilinears can actually generate effective baryon violating trilinears when com- 
bined with the Yukawas (see e.g. [24, 58]). This can be easily seen working at the SO{W) 
level. Take for instance the terms 

W D y/^16^16i,10H + y/^16^10],16H + /i^^lO^lO// +^^10^T6h16h , (B.27) 

where fl'^ = ^/i^Q -|- ??J^q (45//) -|- A^q (54//)^ and by integrating out the pairs IOh — Wh and 
16h — 16// one gets^^ 

hoI^ie^iQUo'^ + h^iQ^iouo' . (b.28) 

/^lO /"16 

After projecting these operators on the MSSM fields, only the second one gives a low-energy 
contribution, leading to the trilinear operator (A^)*"''^ n^d^d^ with 



^X'^y^^ = —Yl'^d^ . (B.29) 



^16 

The SO{W) trilinear operator A*-'''^16^16j;,10^ has no projection on the light MSSM fields 
as well and thus does not contribute to the effective low-energy RPV superpotential. On 



^The argument should be formally carried on at the SM level by integrating out the heavy triplets. 
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the other hand RPV triUnear couplings can arise at the non-renormalizable level from the 
following operator 



Mp 



{liljel + 2 d'iijqk + djd'jul) . (B.30) 



Notice that due to the antisymmetry of the Wiq^ in the SU{5) space the interactions 
in Eq. (B.30) are antisymmetric in the first two generation indices: A^;^^ = — A-j^^. 
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